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Abstract 

' The Debye screening mass for a quark-gluon plasma at high temperature is calcu- 
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. lated to next-to-leading order in the QCD coupling constant from the correlator of two 

00 ' Polyakov loops. The result agrees with the screening mass defined by the location of 

O '. 

, the pole in the gluon propagator as calculated by Rebhan. It is logarithmically sensi- 



Ph! fields. 
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five to nonperturbative effects associated with the screening of static chromomagnetic 



One of the fundamental properties of a plasma is the Debye screening mass mn, 
whose inverse is the screening length for electric fields in the plasma. The definition of the 
Debye mass is conventionally given in terms of the small-momentum {k 0) limit of the 
static {uj = 0) Coulomb propag ator + noo(0, k)]'^ = [Pe(0, k)]'^, where noo(c^, k) is the 
longitudinal photon self-energy function and e{u,k) is the corresponding dielectric function: 



The Debye mass can alternatively be defined in terms of the location of the pole in the static 
propagator for complex k: 



At leading order in the coupling constant, the static self-energy has the simple form noo(0, k) = 
m\) independent of fc, and the definitions (|I]) and (Q) are equivalent. Beyond leading order in 
the coupling constant, they need no longer be equivalent. A fundamental question of plasma 
physics is then this: what is the correct general definition of the Debye mass? 

This question applies equally well to a quark-gluon plasma, where the Debye mass 
describes the screening of chromoelectric fields. Because the coupling constant of quantum 
chromodynamics (QCD) is relatively large, higher order corrections to the Debye mass are 
probably not negligible. The question of the correct definition therefore becomes one of prac- 
tical importance. In the case of QCD, the longitudinal gluon self-energy function Hqq^uj, k) 
is gauge dependent. If the Debye mass is relevant to the screening of chromoelectric fields, 
then it must be gauge invariant. Thus, gauge invariance can be used as a guide to the cor- 
rect definition of the Debye mass. Formal arguments due to Kobes, Kunstatter, and Rebhan 
indicate that, in spite of the gauge- dependence of the gluon propagator, the locations of 
its poles are gauge invariant [|T|. This suggests that (Q) is the correct general definition. 
However, it is important to back up these formal arguments with explicit calculations. 

There have been a number of calculations of the next-to-leading order correction to 
the Debye mass using the conventional definition (|I]) ^. The results are infrared finite but 
gauge dependent. The naive application of the pole definition gives a gauge-dependent 
result, which is also infrared- divergent. However, it was recently shown by Rebhan that, 
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m\) = lim Hoo(0, k). 



(1) 



k^ + Hoo(0, k) = at k'^ 




(2) 



with a careful treatment of infrared effects, the pole definition does in fact give a gauge- 
invariant result. 

Since the calculation of Rebhan involves subtle interchanges of limits, it is desirable 
to have an independent verification of this result. It is also desirable to calculate the next- 
to-leading order Debye mass directly from a gauge-invariant quantity, to provide assurance 
that it is relevant to physical quantities. The simplest such quantity is the correlation 
function of two Polyakov loops. A previous calculation of this correlator to next-to-leading 
nontrivial order by Nadkarni 0] gave results that seemed to be incompatible with simple 
Debye screening. In this Letter, we reexamine this calculation and show that, by careful 
treatment of infrared effects, it can be used to extract the Debye mass to next-to-leading 
order. The result agrees with that obtained by Rebhan from the pole in the gluon propagator. 
These results, together with the general arguments of Ref . |1| , provide compelling evidence 
that the pole definition (0) is the correct definition of the Debye mass beyond leading order 
in the coupling constant. 

We begin by reviewing the calculation by Rebhan. At leading order in the QCD 
coupling constant g, the Debye mass is given by a straightforward perturbative calculation: 



= y C^Nc + Nf)/6 gT, where Nc = 3, Nf is the number of flavors of light quarks, and T 
is the temperature. At next-to-leading order, it is necessary to resum perturbation theory 
by including the Debye mass in the Coulomb propagator. The contribution of order g^ to 
the static longitudinal gluon self-energy is 



5noo(0,fc) = N,g^T J 
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(P + mg,)(p2 + 2p-k) 
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where q = p + k. In (j^), ultraviolet divergences are understood to be removed using 
dimensional regularization. Integrating over p, we obtain the result 
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Upon taking the limit k irriD, we find a correction that is logarithmically divergent and 
gauge dependent. However, as pointed out by Rebhan, for k near imo, the integral in (||) 
is extremely sensitive to the infrared region p — 0. He regularized the infrared region by 



replacing 1/p^ by l/(p^ + m^^^g). Taking the limit k irriD in the presence of the regulator 
and then taking the limit mmag ^ 0, he obtained the gauge- invariant result 
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The logarithm indicates that the Debye mass at next-to-leading order is weakly sensitive to 
the screening of the static chromomagnetic interaction, which is dominated by nonpertur- 
bative effects. The particular regulator used by Rebhan corresponds to a model for these 
nonperturbative effects in which they simply shift the pole in the transverse gluon propaga- 
tor from p = to the imaginary values p = ±zmmag- One could equally well use dimensional 
regularization as the infrared regulator. Again the pole definition gives a gauge-invariant 
result which diverges as the regulator is removed. 

We next review Nadkarni's calculation of the next-to-leading order correction to the 
correlator of Polyakov loops. The calculation was carried out in static gauge, where the 
Polyakov loop reduces to a simple exponential: Q{x) = exp{—igA4{'x)/T). The correlator of 
two such operators separated by a distance R is Cpl{R) = (Trfi(0)Trn(R))/A^^. Nadkarni 
calculated this correlator to order and found 



Cpl{R) 



1 + 



V{Rf 
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where V{R) = e~™'^^/(47ri?) is the Debye-screened Coulomb potential. The complete ex- 
pressions for the functions fi{R) and f2{R) are given in Ref. 0]. As in the Appendix of 
Ref. the integrals can be simplified by expanding out the numerators and cancelling 
propagator factors wherever possible. Keeping only those terms that grow at least linearly 
with R, the functions reduce to 



MR) = -2mnV{Rr'J 
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(27r)3 J (27r)3 (pf + ml){ql + ml){pl + ml){ql + ml){pi - ps) 



where = Pj+k. Ultraviolet divergences are understood to be regularized using dimensional 
regularization. Evaluating the integrals in the limit R ^ oo, Nadkarni obtained 

MR) !^ [-2 log(2m^i?) + 3 - 27] , (9) 
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MR) ^[\ogimnR)+l], (10) 
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where 7 is Euler's constant. The terms proportional to R can be absorbed into an order- (7^ 
correction to the Debye mass in the g'^ term in (j^), which is proportional to exp{—2m£,R). 
The terms proportional to R log R, however, cannot be absorbed into the Debye mass. Nad- 
karni concluded that the next-to-leading order correction to the Debye mass is undefined. 

We now show how the next-to-leading order Debye mass can be extracted from the 
calculation of the correlator of Polyakov loops. The R log R terms in (P) and indicate 



that the integrals in and (§) are sensitive to the infrared region. As a regularization of 
the infrared region, we replace the massless propagator 1/p^ in (|^) by l/(p^ + '^mag); ^^'^ 
similarly for l/(pi — P2)^ in (H)- If we first set mmag = and then evaluate the integrals in 
the limit R 00, we recover Nadkarni's results and (plQl). However if we take the limit 
— > 00 in the presence of the regulator, we obtain a different result. 

We first consider fi{R) in (^. Evaluating the first integral over p with dimensional 
regularization, we obtain the constant value —mo /An. The second integral over p gives a 
function of k that has branch points ai k = ±i{mD + mmag)- The integral over k can be 
evaluated by expressing it as an integral over k along the entire real axis, and then deforming 
the contour into the upper-half complex plane. It receives contributions from the double pole 
at A; = im£, and from the branch cut beginning at A; = i(m£) + mmag)- In the limit R 00, 
the contribution from the branch cut is exponentially suppressed by a factor of e"™™''^^^ and 
can therefore be neglected. The second integral over p in (|^ can therefore be replaced by 
its value at k = irriD, and it reduces to log(2mD/minag)/(8vrmi3) in the limit rrimag ^ 
Evaluating the remaining contour integral over k, we find that the asymptotic behavior of 
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We next consider /2(-R) in (H). In the limit mmag — ^ 0, the integral over pi and p2 
gives log((/c^ + 4m|,)/4mf))/[167r^A;^(A;^ + 4m|,)], which has coincident poles and logarithmic 
branch points aX k = ±2im£,. In the presence of the regulator, the integral has branch 
points at A; = ±2im£, and k = ±i{2mo + mmag)- The integral over k can again be evaluated 
by expressing it as an integral over k along a contour that is deformed into the upper-half 
complex plane. The contribution from the branch cut beginning at = 1(27710 + '^mag) is 
suppressed by a factor of e"'"™''^^ and can be neglected in the limit i? — > 00. The contribution 
from the branch point beginning at k = 2imD has the asymptotic behavior 
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h{R) — ^^logi?. (12) 
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There is a quadratic infrared divergence that is cut off by the magnetic mass, but there are 
no terms linear in R. 



Inserting (|TT]) into (^), we find that the correction proportional to R has precisely 
the form —26m£)R, where ^m^) is given by Rebhan's expression (|^). This correction can 
be absorbed into the order-gf^ term in by making the substitution uio rriD + 6m£, in 
the potential V{R). It therefore represents a correction to the Debye mass. This result was 
obtained by using a magnetic mass as an infrared regulator, which separated the poles asso- 
ciated with the Debye mass from the branch cuts associated with the emission of transverse 
gluons. We could equally well have used a different infrared regulator, such as dimensional 
regularization. As long as the same regulator is used in the calculation of the correlator of 
Polyakov loops and in the calculation of the pole in the gluon propagator, the next-to-leading 
order correction to the Debye mass will be the same. This is evident from comparing the 
relevant parts of the integrals over p in and (0). 

The fact that the Debye mass depends on an infrared cutoff at next-to-leading order 
should not be too disturbing. At large R, the correlator of Polyakov loops is dominated 
by a sum of terms that fall exponentially like exp(— Mi?), with M being the mass of a 



color-singlet bound state in (2+l)-dimensional extended QCD ||5|. This field theory is an 
SU{3) gauge theory describing magnetic gluons, coupled to an adjoint scalar field describing 
electric gluons. One of the bound states consists predominantly of a pair of electric gluons, 
and its mass is 2m/) plus a small binding energy. It is the contribution of this bound state 
that corresponds to the perturbative contribution to the Polyakov loop correlator calculated 
above. The infrared divergence in the next-to- leading order correction to 2m/) is a signal 
that the binding energy of this state cannot be calculated in perturbation theory. 

As an aside, it should be pointed out the truly asymptotic behavior of the Polyakov 
loop correlator as i? — oo is dominated not by the lowest bound state of electric gluons, 
but by a magnetic glueball 0. Its coupling to the Polyakov loop is suppressed by extra 
powers of the coupling constant, but its mass is much less than 2m£) at weak coupling, so it 
dominates at sufficiently large R. 

We have calculated the Debye mass for a quark-gluon plasma to next-to-leading order 
in the QCD coupling constant from the gauge-invariant correlation function of Polyakov 
loops. The result agrees with that obtained by Rebhan from the pole in the gluon propagator. 
It is gauge invariant, but depends logarithmically on nonperturbative effects associated with 
the screening of static chromomagnetic fields. This result provides support for the general 
definition of the Debye mass in terms of the location of the pole in the propagator of the 
gauge field. 
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